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Gamma-ray bursts (GRBs) are some of the most extreme events in the Universe. As well as pro¬ 
viding a natural laboratory for investigating fundamental physical processes, they might trace the 
cosmic star formation rate up to extreme redshifts and probe the composition of the intergalactic 
medium over most of the Universe’s history. Radio observations of GRBs play a key part in deter¬ 
mining their physical properties, but currently they are largely limited to follow-up observations 
of y-ray-detected objects. 

The SKA will significantly increase our ability to study GRB afterglows, following up several 
hundred objects in the high frequency bands already in the “early science” implementation of the 
telescope. SKA1-MID Bands 4 (4 GHz) and 5 (9.2 GHz) will be particularly suited to the detec¬ 
tion of these transient phenomena. The SKA will trace the peak of the emission, sampling the 
thick-to-thin transition of the evolving spectrum, and follow-up the afterglow down to the time the 
ejecta slow down to non-relativistic speeds. The full SKA will be able to observe the afterglows 
across the non-relativistic transition, for ~ 25% of the whole GRB population. This will allow us 
to get a significant insight into the true energy budget of GRBs, probe their surrounding density 
profile, and the shock microphysics. The SKA will also be able to routinely detect the elusive 
“orphan afterglow” emission, from the population of GRBs whose jets are not pointed towards 
the Earth. We expect that a deep all-sky survey such as SKA1-SUR will see around 300 orphan 
afterglows every week. We predict these detection to be >1000 when the full SKA telescope will 
be operational. 
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1. Introduction 


Gamma-ray bursts (GRBs) are among the most extreme events in the Universe, they provide 
a laboratory for investigating fundamental physical processes, trace the cosmic star formation up 
to very high redshifts, and probe the composition of the intergalactic medium over most of the 
Universe’s history. Radio observations of GRBs play a key part in determining their physical 
properties by complementing the broadband spectrum and providing information unique to the 
radio regime. 

The origin of GRBs remained elusive for three decades after their initial discovery (Klebe- 
sadel et al. 1973), but this changed with the discovery of their afterglows at X-ray (Costa et al. 
1997), optical (van Paradijs et al. 1997), and radio (Frail et al. 1997) frequencies. GRBs occur at 
cosmological distances (Meegan et al. 1992; Metzger et al. 1997), and are generally thought to be 
caused by either the collapse of a massive star (Woosley 1993), or the binary merger of two neutron 
stars or a neutron star and a black hole (Eichler et al. 1989; Narayan et al. 1992). The former are 
the progenitors of so-called long GRBs, while the latter are associated with short GRBs. Long and 
short refer to the GRB duration in gamma-rays, and the divide between the two classes is at ~ 2 s 
(Kouveliotou et al. 1993), albeit it should be noted that the measured duration is energy-dependent. 
In both cases, the rapid accretion of material onto the compact object that is formed in these catas¬ 
trophic events is thought to produce a pair of highly collimated relativistic jets, which decelerate 
and eventually evolve into a roughly isotopic blast wave. The isotopic equivalent kinetic energy 
released in a GRB is of order 10 53 ~ 54 ergs, which is similar to that seen in supernova explosions 
when corrected for collimation and radiative efficiency (van Paradijs et al. 2000). However, GRB 
jets have bulk Lorentz factors of order of hundreds, and therefore relativistic beaming makes them 
some of the most luminous events in the Universe. 

By studying the broadband emission of GRBs, from low-frequency radio waves to high-energy 
gamma-rays, we can determine the physical properties of the jet and its surroundings, and also the 
nricrophysics of the emission we observe. GRBs are amongst the highest redshiftobject known: the 
current highest spectroscopic redshift is z = 8.23 (Tanvir et al. 2009), and the photometric record- 
holder is z = 9.2 (Cucchiara et al. 2011). Therefore they can be used as probes of the early Universe, 
of the host galaxies (see e.g. Stanway et al. 2014; Perley & Perley 2013), and of the medium in 
between their host galaxies and us. Furthermore, GRBs are excellent laboratories for studying 
particle acceleration physics, because they are not only accelerating electrons to extreme relativistic 
velocities to produce the observed emission, but are also prime candidates for producing ultra-high 
energy cosmic rays. In this chapter we first focus on the role of radio observations in GRB studies 
(Section |2|) and the progress that has been made in this field over the decade since the original SKA 
science case was written (Section |j). We have performed simulations to give prospects for GRB 
studies in the SKA era, for GRBs whose jets point towards the Earth (Section [4), and also those 
which are pointed away and from which we can detect late-time emission when the GRB outflow 
enters our line of sight (Section [5). For a more comprehensive review on radio observations of 
GRBs and their role in GRB jet studies we refer to Granot & van der Horst (2014). 


2 



GRB SKA 


Davide Burlon 


2. Radio observations of GRBs 

Almost all GRBs have been detected first with 7 -ray instruments and then followed up at other 
wavelengths (see e.g. the interesting case reported by Cenko et al. 2013). Since the discovery of 
the first radio afterglow (Frail et al. 1997) roughly one third of all GRBs with accurate locations 
has been detected at radio frequencies (Chandra & Frail 2012). This rate is much lower than at 
higher observing frequencies, where for instance 93% of GRBs detected in gamma-rays by the 
Swift satellite are also detected at X-ray frequencies and 75% have been detected in optical. This 
difference is due to the fact that GRBs have radio flux densities typically at the sub-mJy level, 
which for a long time has been close to the sensitivity limits of even the largest radio observatories. 
Besides the sample of radio detected GRBs being sensitivity limited (Chandra & Frail 2012), it has 
been claimed that there are in fact two populations of radio-bright and radio-faint GRBs (Hancock 
et al. 2013). The latter has been found by stacking radio visibility data of many GRB follow-up 
observations. This could be tested with the upgraded Very Large Array and ultimately with the 
SKA. If there are indeed two populations of radio afterglows, the cause of this can be resolved by 
broadband modelling of a large sample of individual GRBs. 

In the standard model for GRBs, the afterglow is thought to be emission from the blast wave at 
the front of the relativistic jet pointing towards us. As this shock sweeps up the ambient medium, it 
accelerates electrons (and protons) and amplifies the magnetic field, both necessary for the observed 
synchrotron radiation. Deceleration of the blast wave results in the peak of the spectrum moving 
from high to low observing frequencies over time (Sari et al. 1998). This spectral peak is typically 
already below optical frequencies when the first observations commence, resulting in declining 
light curves at optical and X-ray frequencies, but rising light curves in the radio at early times, 
up to the peak occurring days or weeks after the initial gamma-ray trigger (see Figure [lj). Multi¬ 
frequency radio observations at later times can follow the evolution of both the peak frequency and 
peak flux. Furthermore, the radio regime is affected by synchrotron self-absorption, which makes 
early radio detections challenging, but when detected, time evolution of the self-absorption fre¬ 
quency provides important constraints on the blast wave physics. Therefore, the radio regime plays 
an important role in constructing the full broadband spectrum, to constrain both the macrophysics 
of the jet, i.e. the energetics and ambient medium density, and the microphysics of the electrons 
and magnetic fields necessary for synchrotron emission (Wijers & Galama 1999). 

Radio afterglow emission peaks much later than at higher frequencies, and it also lasts much 
longer, for months or even years (Figure [I]; e.g., Frail et al. 2000; Berger et al. 2004; van der 
Horst et al. 2008). Late-time observations can be used to study the GRB jet as it decelerates into 
the non-relativistic phase. This provides a unique probe for the dynamics of the jet, in particular 
for determining the true total kinetic energy in the blast wave; at late times the outflow is almost 
spherical and energy estimates are not complicated by relativistic effects (see Section 0 for a more 
detailed discussion). In addition to these late-time studies, the radio regime is also well suited 
for studying early-time phenomena which may go undetected at high frequencies, for instance 
emission from the reverse shock. When the blast wave at the front of the jet is formed, there is 
not only a shock moving forward into the ambient medium, but also a shock moving backwards 
into the ejecta (Sari & Piran 1999). Under the right conditions this reverse shock will result in an 
optical flash in the first tens of seconds after the gamma-ray trigger (Akerlof et al. 1999), which can 
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GRB970508 GRB980703 GRB030329 



Time (days after burst) 

Figure 1: Radio light curves at 4.9 and 8.5 GHz (top panels) and spectral indices (bottom panels) for the 
well-sampled, long-lasting GRBs 970508, 980703 and 030329 (Granot & van der Horst 2014). The spectral 
index a (where F v °e v“) between 4.9 and 8.5 GHz varies significantly due to the spectral evolution and 
scintillation effects. The dashed lines in the bottom panels indicate spectral indices of 2, 1/3 and —0.6, 
which are expected below the self-absorption frequency, between the self-absorption and peak frequency, 
and above the peak frequency, respectively. 


only be detected if robotic telescopes are on target fast enough. However the peak of this emission 
moves to lower frequencies over time and can be probed at radio frequencies on a time-scale of 
hours to days (Kulkarni et al. 1999). Multi-frequency observations of this reverse shock provide 
crucial information regarding the particle content of the jet and its internal magnetic field (e.g. 
Perley et al. 2014; Anderson et al. 2014, and references therein). 

Another way in which radio observations play a key role in GRB studies is by providing 
evidence for the relativistic expansion of the jet. The source size and its evolution can be estimated 
either directly by means of Very Long Baseline Interferometry (VLBI), or indirectly by utilising 
the effects of interstellar scintillation. The former method has only been possible for one nearby 
GRB (Taylor et al. 2004), and the source size measurements have been combined with its light 
curves to better constrain the physical parameters (Granot et al. 2005; Mesler & Pihlstrom 2013). 
Indirect source size measurements are possible because of scintillation due to the local interstellar 
medium, which modulates the radio flux (Goodman 1997). While the angular size of the blast 
wave is initially smaller than the characteristic angular scale for interstellar scintillation, it grows 
with time until it eventually exceeds this scale and the modulations quench. This can be utilised in 
determining the source size and the expansion speed of the blast wave, which impose constraints on 
the modelling of GRB afterglows (Frail et al. 1997). The future possibilities of VLBI observations 
of GRBs in the SKA era arc discussed elsewhere in this book; see the Chapter by Paragi et al., 
“Very Long Baseline Interferometry with the SKA”, in proceedings of “Advancing Astrophysics 
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with the Square Kilometre Array”, PoS(AASKA14)143. 

All the observational characteristics at radio frequencies discussed so far are for GRBs in 
which the jet is oriented along our line of sight. Since the initial Lorentz factors of the shocks arc of 
the order of a few hundred, relativistic beaming causes the initial y-ray emission to go undetected 
if we arc not observing a GRB within the jet opening angle. Therefore, we only detect a small 
fraction of the total number of events (Rhoads 1997). Estimates of collimation for individual GRBs 
arc usually derived from the afterglow light curve: an achromatic break to a steeper decay indicates 
a transition across the jet opening angle. The time at which this ‘jet break’ occurs, combined with 
the ambient medium density and the energetics of the GRB, leads to an estimate of the jet opening 
angle (Sari et al. 1999). This approach is limited by the difficulty in identifying an unambiguous 
jet break in the light curve and yields a wide variation of jet angles, but typically a few to a few 
tens of degrees. An alternative way to obtain a handle on jet opening angles, and on the total GRB 
rate whether beamed towards us or not, is to carry out a survey for GRB afterglows at very late 
times when the outflow is (quasi-)spherical and only detectable at radio frequencies. Since a large 
fraction of GRBs arc associated with supernovae of Type Ib/c, it has been suggested to search for 
radio counterparts of relatively nearby supernovae of that type (Paczynski 2001; Granot & Loeb 
2003). Attempts to find such orphan radio afterglows have been unsuccessful thus far (Berger et al. 
2003; Soderberg et al. 2006a; Bietenholz et al. 2014). 


3. Progress in the Last Decade 

Our understanding of GRBs has increased significantly over the past decade, mainly because 
of improvements in observational capabilities across the electromagnetic spectrum. For instance, 
the Swift satellite has revealed very interesting X-ray behavior in the first minutes to hours (Nousek 
et al. 2006), robotic telescopes have done the same at optical frequencies (e.g., Racusin et al. 
2008), and the Fermi Gamma-ray Space Telescope has opened the high-energy gamma-ray regime 
for GRBs, including emission at GeV energies (Ackermann et al. 2013). Furthermore, the GRB- 
supernova association has been further established for long GRBs (Woosley & Bloom 2006), while 
afterglows of short GRBs have been discovered (Gehrels et al. 2005), and strong indications for 
their connection with binary mergers of compact objects have been found (Tanvir et al. 2013). 

From the radio perspective, GRB 030329 had a major impact: the first high-luminosity GRB 
at low redshift, which also satisfied the GRB-supernova connection. It is the only GRB up to now 
for which the source size has been measured with VLBI, and the radio afterglow was so bright and 
long-lasting that it has been detected for almost a decade at low radio frequencies (van der Horst 
et al. 2008; Mesler et al. 2012). The multi-frequency light curves have made it possible to follow 
the evolution of the blast wave from the ultra-relativistic to the non-relativistic phase and perform 
broadband modeling in the different phases. This has led to tight constraints on the physical pa¬ 
rameters, and also a consistent picture together with modeling of the VLBI measurements (Granot 
et al. 2005; Mesler & Pihlstrom 2013). The puzzle that still remains is that all the modeling work 
points to the ambient medium being homogeneous, and in that case a signature from the counter-jet 
(the jet originally pointed away from the Earth) should have been observed (De Colle et al. 2012). 
This did not happen, questioning the existance and properties of this counter-jet. 
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Figure 2: The left panel (Anderson et al. 2014) shows the radio light curves of the bright GRB 130427A 
at 5 and 15 GHz (see also Laskar et al. 2013; Perley et al. 2014; van der Horst et al. 2014), showing the 
evolution of the radio peak moving from high to low frequencies at early times. The right panel shows 
strong variability at the peak of the 5 GHz light curve on a timescale of hours (solid symbols), and perhaps 
even shorter (15 minutes; open symbols), due to interstellar scintillation (van der Horst et al. 2014). 


Radio observations have played an important role in the broadband modeling of a larger sample 
of GRBs, e.g., very energetic GRBs detected by Swift (Cenko et al. 2010) and those for which 
high-energy gamma-rays have been detected by Fermi (Cenko et al. 2011). These studies indicate 
that there is a large spread in the physical parameters, which has implications for models of GRB 
progenitors, particle acceleration, and magnetic field amplification by shocks (Granot & van der 
Horst 2014). A striking example is GRB 130427A, another nearby, high-luminosity event with 
extremely good temporal and spectral coverage, over 10 orders of magnitude in time and 16 orders 
of magnitude in observing frequency (Ackermann et al. 2014; Maselli et al. 2014). A bright radio 
flare was observed in the first days, with the peak moving from high to low radio frequencies 
over time (left hand panel of Figure [2j; Laskar et al. 2013; Anderson et al. 2014; Perley et al. 
2014; van der Horst et al. 2014). The earliest clear detection at 8 h and subsequent observations 
at high cadence have shown how early-time radio observations can contribute significantly to our 
understanding of the physics of both the forward and reverse shocks. This was also the second 
GRB for which radio variability on the timescale of hours has been reported (right hand panel 
of Figure |j; van der Horst et al. 2014), and in this case, as well as for GRB 070125 (Chandra 
et al. 2008), the variability strength and timescale are both consistent with interstellar scintillation 
effects. 

An important probe for getting a better understanding of the jet and magnetic field structure 
is polarization measurements. In the optical, linear polarization of the order of a few percent has 
been observed in emission from the forward shock (Covino et al. 1999; Wijers et al. 1999), while 
this can be tens of percent for the reverse shock (Mundell et al. 2013). Optical circular polarization 
has recently been discovered, at a level of only a few tenths of a percent (Wiersema et al. 2014). At 
radio frequencies the most stringent constraints have been obtained for GRB 030329 at late times 
(Taylor et al. 2005), and for GRB 130427A at early times (van der Horst et al. 2014). The upper 
limits are only a few percent in both linear and circular polarization, and thus getting close to the 
optical detections. Deeper radio polarization measurements for bright GRBs, preferably at early 
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times and for GRBs with optical polarization detections as well, will provide more insight into the 
magnetic field strength and structure inside the GRB jet. 

Coherent radio emission has been predicted at earlier times than from the reverse shock, which 
peaks at a timescale of hours to days at radio frequencies. In these models the jet is magnetically 
dominated (Usov & Katz 2000; Sagiv & Waxman 2002; Moortgat & Kuijpers 2006), and detecting 
this emission would uniquely constrain the jet production and collimation mechanism. The first 
searches have not resulted in any firm detections (Bannister et al. 2012; Obenberger et al. 2014), 
but surveys with large field-of-view telescopes are vital for this research. Such studies should also 
clarify if (some of) the recently discovered fast radio bursts (Lorimer et al. 2007; Thornton et al. 
2013) are connected to GRBs (Zhang 2014). 

Radio surveys with sufficient sensitivity and sky coverage are not only important for prompt 
radio emission searches, but also for detecting orphan afterglows. Such studies would let us directly 
constrain the overall beaming fraction of GRBs, and yield the rate of the events that generate 
this phenomenon. Furthermore, since radio emission is unaffected by extinction, and a significant 
fraction of GRBs arc thought to occur in highly obscured regions (Perley et al. 2009), radio surveys 
can provide a new picture of the GRB population unbiased by their local environments. This ties 
in with studies of dark bursts, i.e. GRBs for which the optical afterglow is much dimmer than 
expected from the observed X-ray emission. These GRBs arc either at a high redshift or there is 
significant optical extinction in the host galaxy. The latter is the most common explanation, and the 
amount of extinction can be estimated by modeling of the X-ray and radio light curves (Rol et al. 
2007; Zauderer et al. 2013). 

4. The GRB population pointing to Earth 

In this section we focus on GRBs that can be studied by the SKA after a trigger has been 
produced by y-ray detectors, as now happens with Swift and Fermi providing the high energy 
trigger and initiating the follow-up campaign. These arc the GRBs that have their jets pointed 
towards the Earth. 

Substantial advances in breaking the degeneracies of the physical parameters of the afterglow 
emission process and, globally, of the GRB itself can only be possible through a considerable 
increase of the detection rate of GRBs coupled to the ability to follow the radio emission from 
early to very late times (i.e. at extremely faint flux densities). 

The key questions we will address arc: 

1. What are the prospects for the detection and follow up of GRBs by the SKA? 

2. What is the expected rate of GRBs that can be observed by the SKA? 

3. What is the best observational strategy, considering the wide frequency range covered by the 
SKA and the behaviour of the GRB emission at radio frequencies as a function of time? 

In order to answer these questions, we first characterised the radio afterglow light curves of a 
synthetic population of GRBs constructed so that its y-ray properties (flux distributions, energet¬ 
ics/luminosities) match those of the real populations of GRBs detected by CGRO/BATSE, Swift 
and Fermi. GRBs are distributed in redshift between z = 0 and 10 according to the GRB formation 
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rate (see Salvaterra et al. 2012, for the different components that are used). For each simulated 
burst, the radio afterglow emission was simulated using a hydrodynamic numerical code (for full 
details we refer the reader to van Eerten & MacFadyen 2011, and van Eerten et al. 2012). The 
micro-physical parameters (how much energy is distributed to the elections and magnetic fields, 
the slope of the electron distribution, and the density of the homogeneous ambient medium) were 
set in order to reproduce the typical fluxes of the afterglows observed in the optical (Melandri 
et al. 2014) and radio (Ghirlanda et al. 2013a) bands of a carefully selected sample of bright GRBs 
observed by Swift (the so-called BAT6 sample of Salvaterra et al. 2012). These are the brightest 
GRBs detected by Swift/BKT and they have typical flux densities between 1 and 5 days distributed 
between 0.1 and 2 mJy at 8.4 GHz (Chandra & Frail 2012). The micro-physical parameters were 
tuned to reproduce these flux levels (as shown by the red points in the right panel of Fig. [?). This 
approach is discussed in detail by Ghirlanda et al. (2013b). A caveat is the lack of the possible 
contribution of the reverse shock component, which might dominate the emission at early times. 
Estimates on the detection rates at early times should hence be conservative lower limits. Another 
caveat is that the simulations do not consider the possible supernova contribution, nor the effects 
of scintillation. 
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Figure 3: Left: Cumulative flux distribution of afterglows at SKA1-LOW and SKA1-MID observing bands 
2, 4, and 5 (from left to right). The shaded area represents the span between 2 and 10 days (left and right 
boundary respectively). The dashed line represent the flux at the peak time. Right: Afterglow peak flux vs. 
observed time of the peak (which happens at the time the peak of the spectrum crosses the observed band) for 
different SKA frequencies. The colours (red, blue, black, and green) represent bands 5, 4, 2 of SKA1-MID 
and SKA1-LOW, i.e. 9.2, 4, 1.4 and 0.2 GHz respectively. 

At early times (days after the burst), the emission from the forward shock is suppressed by 
synchrotron self-absorption, and therefore the flux density due to this component is expected to 
rise until a spectral peak is produced by the passage of the frequency (be it the injection frequency 
or the self-absorption frequency) corresponding to the peak of the spectrum in the telescope ob¬ 
serving band. This effectively means that the observing strategy should consider the available 
frequencies and the typical behaviour of the flux of the afterglow at those frequencies. For in¬ 
stance, the detection rate is maximised if the observations are performed close to the time when 
the peak of the afterglow is reached at a certain frequency. On the other hand, early time observa¬ 
tions, which are critical for the multi-wavelength study of the GRB emission (coupling radio with 
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optical/NIR, and X-rays), should be planned considering that the radio spectrum of the afterglow 
is still self-absorbed. Therefore, a successful detection is optimised if performed at relatively high 
radio frequencies 1 . 

In Fig. |] (left) we show the cumulative flux distribution of the simulated sample of GRBs. The 
different colours represent different bands of the SKA. The shaded regions represent the cumulative 
flux distributions if observations are performed between 2 and 10 days after the GRB trigger. The 
peak flux distributions are shown by the dashed lines. Two behaviours are clearly visible. Firstly, 
the highest flux is always at the time of the second observation (i.e. the emission is still in the 
self-absorbed region of the spectrum). Secondly, the width of the region decreases with increasing 
frequency, because observing in a higher radio band will bring the observation closer to the region 
in the spectrum where most of the flux is emitted, and where the evolution is faster. SKA 1-MID 
band 5 (9.2 GHz) is clearly the band that tends to pick a radio afterglow between 2 and 10 days 
closest to the time of the peak. Another way of looking at the distance of the dashed line from 
the corresponding shaded area (which again, is maximum for SKA 1-LOW, and minimum for the 
highest SKA1-MID band shown), is a measure of how long to delay an observation in order to 
maximise the chance of detection in each SKA band. 

Table 1: SKA detection rates of radio afterglows independent from high-energy triggers. For each band we 
report the 5 <7 sensitivity limit. A standard observation of 12 hours is assumed. " Refers to “early science”, 
an initial implementation with half the sensitivity across all bands. b The sensitivity is assumed to be 4 x 
better in the LOW band, and 10 x otherwise. 




SKA1-MID 


SKA 1-MID E.S." 

SKA2 fc 

Band 

Freq. 

[GHz] 

Flim 

[MJy] 

Rate 

[# sr~ 1 yr" 1 ] 

Rate 

(# sr -1 yr^ 1 ] 

Rate 

[# sr- 1 yr" 1 ] 

LOW 

0.2 

3 

0 

0 

2 

Band 2 

1.4 

1.39 

300 

170 

465 

Band 4 

4.0 

0.65 

465 

450 

465 

Band 5 

9.2 

0.75 

465 

465 

465 


Alternatively, one can visualise this behaviour by investigating the flux density at which a GRB 
afterglow peaks when the transition from optically thick-to-thin happens. Fig. |] (right) shows a plot 
of the peak fluxes versus the time at which the frequency corresponding to the peak of the spectrum 
crosses the band (different colours represent the simulated SKA bands). The more energetic GRBs 
peak at later times, with the lower frequencies lagging the high ones. The brighter bursts populate 
the upper right part of the distributions, and can indeed be brighter than 0.1 mJy at 1.4 GHz even 
months after the prompt emission has disappeared. However they are very rare, representing less 
than 1% of the whole population. We note that this is probably a lower limit on the actual fraction 
of these highly luminous afterglows, given that small number statistics dominate the tail of the 
simulated distributions. We have calculated the rates of GRB afterglows that we expect to be able 
to observe at the time of their peak, depending on the level of sensitivity 2 the SKA1 and SKA2 

1 Conversely, the chance of detection of the reverse shock component is maximised in the first day. 

2 Sensitivities are calculated from the SKA1 system baseline design documents, see Dewdney et al. (2013). 
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will reach in each band. Table [I] summarises these results for the whole synthetic population of 
afterglows (i.e. regardless of their y-ray properties). Given that the logN-logF is almost flat at 
sub-jUjy levels, in particular for the high frequency bands, there is no big improvement in expected 
rates from SKA1 to SKA2. The catch is that in order to detect a GRB pointed to Earth, we will 
always need a trigger. On the other hand, for the GRBs that we will be able to detect, the SKA will 
be able to sample to lightcurve down to the /./ .1 y level, producing finely sampled lightcurves, that 
current facilities are able to produce just for the very close-by GRBs. 

In the past, the observation strategy for GRB afterglows was driven more by the goal of having 
coeval multi-wavelength observations across the spectrum. The need of being on target as soon as 
possible in the optical, NIR, and X-rays possibly drove the attempt of having a detection in the 
very first few days. A small number of bursts were detected and followed-up to extremely late 
times (see Frail et al. 1997, 2004; Berger et al. 2004; van der Horst et al. 2005, 2008), but the 
large majority (> 70%) of radio observations resulted in non-detections overall. Note that when 
just the (y-ray) brighter population of bursts was considered, the ratio of positive detections rose 
to ~ 50% (Ghirlanda et al. 2013a), which hints that the SKA will indeed be able to probe the 
hidden population of radio afterglows (Anderson et al., 2014 -in preparation- will present the first 
systematic radio follow-up of Swift GRBs at GHz frequencies). 

4.1 GRB calorimetry: accessing the true energy budget 

One of the key characteristics of the GRB phenomenon that has eluded direct measurements is 
the estimate of the true energy budget of the central engine that powers the explosion. The isotropic 
equivalent energy emitted in y-rays (E)q so ) is just a proxy of the true GRB energy. Indeed, due to 
the collimated structure of GRBs, the energy E y emitted just in y-rays is a factor ~ ^jet - /2 smaller 
than that derived assuming isotropy, where 6>j et is the jet angle. Moreover, the measure of the jet 
angle can be maned by observational biases and possibly by additional components at play in the 
external shock emission. Additionally, the energy emitted in the high-energy band is a fraction of 
the kinetic energy (£kin,iso), initially available. 

According to the standard model, the fireball becomes non-relativistic (NR) at a time 
which is long after the prompt emission has faded away. The theoretical prediction of when this 
transition should occur is subject to debate, and consequently there is roughly an order of magni¬ 
tude in the range of values that ?nr can take, from hundreds (e.g. Livio & Waxman 2000; Sironi & 
Giannios 2013) to thousands of days (e.g. Piran 2004; Zhang & MacFadyen 2009; Wygoda et al. 
2011; van Eerten et al. 2012). For this Chapter, we have adopted the “best case” scenario, with 
the caveat that this is subject to further improvements. Following Fivio & Waxman (2000), the 
transition happens at ?nr r*j 275(1 +z)(n /1 cm- 3 )- | / 3 (£ kin;iso /10 53 erg) 1 / 3 days. The advantage of 
performing energetic estimates (calorimetry) at such late times is that the afterglow emission be¬ 
comes independent from the bulk Forentz factor (which becomes ~ I). This allows us to derive the 
true GRB kinetic energy (e.g. Frail et al. 2000; Shivvers & Berger 2011), which is directly related 
to the energy emitted in y-rays during the prompt phase through a radiative efficiency factor t] . This 
is related to the measurable Eyj so by Ftinjet^ [$jet 2 /(2rj)] x£ 7 j so . This is where the power of these 
observations becomes evident: assuming the efficiency tj of the shocks will give us an independent 
measure of the jet opening angle; conversely for the cases in which there will be a measure of the 
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jet opening angle (e.g. from the achromatic break in the X-ray and optical afterglow), we will be 
able to compute the shock efficiency. 
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Figure 4: Transition to the NR regime: the flux is calculated at time of the crossing to the NR regime, 
calculated as per Livio & Waxmann, 2001. The colour code is as per Fig. 3. 


In addition, observations before and after NfR are probes of the ambient medium profile. In 
the simple case of a homogeneous density ambient medium, a peculiar flattening of the lightcurve 
is predicted by the standard model. In the simple case of uniform density ISM, a flattening of 
Aa = (21 — 5p )/10 is expected for the lightcurve across this transition (Frail et al. 2005), for 
the standard values of the electron distribution slope p. Interestingly, a steepening could also 
potentially be observed at times after ?nr, and that would also help us constrain the microphysical 
parameter p. In the case of a wind-like medium, the transition is expected at later times, with a 
distinct slope variation of the lightcurve. Therefore observing across the NR transition is a powerful 
diagnostic of the circumburst medium density profile. 

We have computed the expected flux density at the NR transition for the SKA bands 2, 4, and 
5, see Fig. |4j. Consistently with the observations available so far, only a handful of objects reach the 
sensitivity limit of current facilities (a few p)y at best). We foresee that SKA1-MID will be able 
to observe a minor fraction of these (of the order of 2-5%, down to 0.9-2% for the “early science” 
phase), but that the full SKA will routinely observe a significant fraction (15-25%) of the whole 
GRB afterglow population at the NR transition.Even for the full SKA though, a substantial (75- 
85%) fraction of sources will be undetectable. We calculated the fractions based on the number of 
afterglows above a 12-hour-long observation (5a), over the total population. 

4.2 The case for a very sensitive high-energy telescope 

The sensitivity of the SKA will allow us to observe almost the complete population of GRBs, 
provided that a y-ray instrument, more sensitive than Swift will be operational in the SKA era, 
and provided that a GRB localisation to a few arcsec spatial resolution will be available. We have 
explored the case of a burst monitor with a (~ 5x) greater sensitivity than that of the Swift/BAT. 
As can be seen in Fig. |5j, the bursts that can be seen in band 5 (9.2 GHz) of SKA 1-MID are a larger 
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fraction of the whole population with respect to the ones that will trigger Swift provided that a 
7 -ray trigger is supplied. 



Figure 5: LogN-logF of SKA1-MID band 5 (9.2 GHz) showing the flux at the peak of the light curve for the 
population of simulated GRBs. The red solid and blue dashed lines represent the population of GRBs that 
have a 7 -ray photon peak flux (P) greater than 0.08 and 0.4 ph cm ~ 2 s' 1 , respectively. The latter represents 
the synthetic population of bursts that can trigger Swift. The improvement of radio flux sensitivity assured by 
the SKA can be exploited only if a 7 -ray detector more sensitive (in the case shown, we assume a five-fold 
improvement in trigger sensitivity) than Swift will be operating in the SKA era. The distance between the red 
and blue line shows how many GRB afterglows, whose radio flux density is accessible by the SKA, we will 
miss out unless there is a GRB trigger. The black line shows all GRBs pointing to the Earth, irrespectively 
of the their 7 -ray flux. 


In short, the future absence of a burst monitor with an improved sensitivity than that of the 
Swift/ BAT, would imply that the rate of GRB afterglows that are observable at the /r-Jy level is 
comparable to the the ones we can already observe today (i.e. ~ 50 sr~ J yr 1 ). In conclusion, while 
trigger-dependent detection rates might not dramatically improve in the SKA era, the sampling of 
the lightcurves of the GRBs that will produce a 7 -ray trigger will be far better than that available 
today. 

5. Revealing the orphan afterglow population 

Most GRBs do not have their jets pointing towards us. As a result they are undetectable in 
the 7 -ray band but should be detectable as ‘orphan afterglow (OA)’. Given the typical jet opening 
angles of GRBs, the population of OAs should outnumber the pointing GRB population by two 
orders of magnitude. 

Jet opening angles have been derived by the afterglow light curve break (Frail et al. 2001; 
Ghirlanda et al. 2004) and have a distribution clustered at 6>j et ~ 5 deg. It is expected that, purely 
by geometric arguments, for each GRB pointing towards the Earth that can potentially trigger a 
7 -ray instrument, there are ~ 2(0j et /O.l ) -2 ~ 260 events directed in all other directions in the sky. 
More realistically, Ghirlanda et al. (2014) have estimated that — given that 0jet has a distribution 
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of values — the ratio of OAs to GRBs pointing to us is ~ 40. These bursts are undetectable in 
the high energy band because are viewed at angles 0 v j ew > 0j et , but are expected to be potentially 
observable at optical and radio wavelengths from the time when, during the afterglow evolution, 
their Lorentz factors slows down to T ~0 v iew _1 - These events should therefore lack their prompt 
emission (and hence have been dubbed ‘orphan afterglows’) and should show up in deep optical 
and radio transient surveys. The characterisation of their properties, the feature of their - lightcurves, 
the typical timescales and peak fluxes are of particular interest for the deep, large surveys of the 
future facilities like SKA1-SUR, since this class of transients could (i) be difficult to disentangle 
from other classes of transients, and (ii) constitute a considerable fraction of transient events in 
the sky. Up until today the quest for then - discovery has resulted in upper (6>j et ,< 22° Levinson 
et al. 2002) or lower (0j et > 0.8° Soderberg et al. 2006b) limits on the GRB opening angles. Some 
peculiar transients in past surveys could have been genuine OAs both in the optical (Rau et al. 
2006; Zou et al. 2007; Malacrino et al. 2007) and radio (for a recent reveiw on blind surveys see 
Muiphy et al. 2013) regime, but none have been confirmed so far. 



Figure 6: Cumulative flux distribution of orphan afterglows at the frequencies of SKA1-MID/SKA1-SUR 
and SKA 1-LOW are shown. To compute the number of expected OA detections, as detailed in the text, we 
assumed a deep (1000 hour), all-sky (3k) survey. The rates can be rescaled easily since they are linearly 
proportional to the sky fraction covered. 


The evolution of the OA lightcurve differs from that of a normal afterglow before reaching 
its peak flux. In fact, an OA should start to be potentially observable 3 when the edge of the jet 
intercepts the line of sight of the observer: this happens when the fireball has decelerated to F = 
1 / sin(0 v j ew — 0j et ). In this respect, the OA is driven by geometry because it is a dynamical peak, not 
a spectral peak generated by the passage of the peak of the spectrum through the observing band. 
This happens at T = 1/sin(0 v j ew ), and the afterglow evolves like that of a standard GRB thereafter. 
Ghirlanda et al. (2014) have estimated that the distribution of peak times is very broad and peaks at 
~ 1000 days. In addition, we show in Fig. || the cumulative flux distribution of a simulated sample 
of OAs in the SKA1 bands. The rates are expressed in units of deg~ 2 yr 1 , so that the conversion 

3 Under the simplified assumption of a sharp edged cone shape for the jet. 
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to the prediction for different fields of view will be straightforward. The high (low) flux tail of the 
logN-logF are consistent with a slope of —1.7 (—0.4) roughly in all bands. The GRBs with small 
opening angles dominate the logN-logF in the high flux tail; conversely the low flux tail is due to 
the cumulative contribution of GRBs with progressively large opening angles. We calculated from 
Fig. |6|how many OAs an all-sky survey with SKA1-SUR should be able to detect each week in 
every PAF band. This result is summarised in Table |2[ 


Table 2: SKA1-SUR detection rates of orphan afterglows in the SKA1-SUR bands. The rates are derived 
from the flux density disttibutions shown in Fig. 6, assuming 1000 hours of integration time, all-sky (3k), 
5(7 and given in units of 37T 1 week -1 . h] lm is in units of liJv. " Refers to “early science”, an initial 
implementation with half the sensitivity. b The sensitivity is assumed to be 10 x better for SKA2. 




SKA1-SUR 


SKA1-SUR E.S." 


SKA2* 


Band 

Freq. [GHz] 

Flim 

Rate 

Flim 

Rate 

Flim 

Rate 

PAF band 1 

0.7 

1.54 

357 

3.08 

178 

0.2 

1190 

PAF band 2 

1.4 

0.92 

387 

1.84 

297 

0.9 

1011 

PAF band 3 

4.0 

1.23 

297 

2.46 

237 

0.6 

892 


The possibility of revealing the population of oiphan afterglows by SKA through its survey is 
tantalizing. Our estimates suggests that at radio frequencies orphan afterglows should be relatively 
long-term transients with typical duration above the expected flux limit of the SKA survey of hun¬ 
dreds of days. The challenge will be to distinguish the population of oiphans from other possible 
sources (non-GRBs) that produce such slow transients. The unique possibility offered by the SKA 
multifrequency observations will be to study the radio spectrum of these transients. Moreover, as 
discussed more in detail in the Chapter on the VLBI capabilities of the SKA (Paragi et ah, “Very 
Long Baseline Interferometry with the SKA”, in proceedings of “Advancing Astrophysics with the 
Square Kilometre Array”, PoS(AASKA14)143), the longest baselines will secure sub-mas reso¬ 
lution. This is turn will mean that, at least for the OAs with sufficient signal-to-noise we will be 
able to measure the expansion velocity of the ejecta. This will be an additional tool to distinguish 
transients powered by a central engine from radio SNe. 


6. Conclusion and the road to the SKA 

The SKA era will be transformational for the study of GRBs, because it will allow us to 
observe close to 100% of the detected GRBs directed towards the Earth. In addition to very high 
detection rates, if compared to the state of the art, it will be able to produce well-sampled radio 
lightcurves, bringing radio-astronomy to the same level of completeness obtained in the last decade 
at higher (optical, X-ray) frequencies. 

For a quarter of all GRBs with a jet pointing to the Earth, the SKA will be able to trace 
the lightcurve up to very late times, even across the non-relativistic transition. We will get an 
unprecedented insight into the true energy budget of GRB engines by means of these very late-time 
observations. Additionally, these observations will probe both the macrophysics of the ambient 
medium (e.g. the density profile) and the microphysics of the shocks, which are believed to be at 
the basis of afterglow emission. 
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The ability to both increase the sheer rate of detections, and to follow the afterglow lightcurve 
to very late times will naturally depend on the availability of a y-ray trigger with arcsecond resolu¬ 
tion. 

Finally we will be able to observe the so far elusive class of “orphan afterglows”, which arc 
predicted to exceed the standard afterglow rate by orders of magnitude. The full SKA telescope 
should detect up to >1000 of these slowly evolving transients every week. 

On the way to the SKA1, the “early science” phase will roughly triple the rate of observed 
afterglows from GRBs pointed towards Earth. The next steps forward will be taken by surveys 
on SKA pathfinder instruments such as the Australian SKA Pathfinder (ASKAP; Johnston et al. 
2007). The VAST survey (Murphy et al. 2013) at GHz frequencies will be first of its kind to reach 
a uniform sub-mJy sensitivity in the whole southern hemisphere. VAST will not have the sensitivity 
to target the bulk of the afterglow population, but it will be the benchmark for the SKA era blind 
surveys of slow transients. 
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